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, Abstract 

We investigate the confinement-Coulomb phase transition in the 
four-dimensional (4D) pure compact U(l) gauge theory on spherical 
^ ■ lattices. The action contains the Wilson coupling (3 and the double 

q_) \ charge coupling 7. The lattice is obtained from the 4D surface of 

the 5D cubic lattice by its radial projection onto a 4D sphere, and 
made homogeneous by means of appropriate weight factors for indi- 
^ \ vidual plaquette contributions to the action. On such lattices the 

' two-state signal, impeding the studies of this theory on toroidal lat- 

tices, is absent for 7 < 0. Furthermore, here a consistent finite-size 
scaling behavior of several bulk observables is found, with the correla- 
tion length exponent v in the range v = 0.35 — 40. These observables 
include Fisher zeros, specific-heat and cumulant extrema as well as 
pseudocritical values of (3 at fixed 7. The most reliable determination 
of v by means of the Fisher zeros gives v = 0.365(8). The phase tran- 
sition at 7 < is thus very probably of 2 nd order and belongs to the 
universality class of a non-Gaussian fixed point. 



1 Introduction 



1.1 Motivation 

The introduction of a space-time lattice for quantum field theories serves sev- 
eral purposes. It provides a regularization for the renormalization scheme; 
it allows us to apply efficient computational methods to perform the func- 
tional integrals; it may be considered a mere approximation scheme for the 
problem in the continuum. Among the four-dimensional (4D) gauge field 
theories with Lie groups the one with £7(1) symmetry at first sight appears 
to be the simplest to test this approach. It is also of considerable interest as 
it has many features in common with QCD, like a confining strong coupling 
phase, topological excitations and gauge-balls. In addition it shows a phase 
transition (PT) to a phase with weak coupling signature, a massless photon 
and long range interaction. In fact it has been the first lattice gauge model 
with continuous gauge group to be studied by the computational methods 
that became available in the 80's [1.. 

Below (in subsect. II. 2j) we will discuss the various results obtained since. 
However, in summary we may say that up to now there is no definite answer 
to the critical properties of its PT. In most simulations a two-state signal 
at the PT indicated a 1 st order transition. On the other hand, the critical 
behavior according to such a transition has not been confirmed in thorough 
finite-size scaling (FSS) studies. The problem persisted when the original 
Wilson action containing only cos(6p) was extended to include the double 
charge coupling 

S = -J2((3cos(Q P ) +7cos(26 P )) . (1) 
p 

Here 9p 6 [0, 2ir) is the plaquette angle, i.e. the argument of the product 
of U(l) link variables around a plaquette P. It was conjectured, that the 1 st 
order PT changes into a 2 nd order one at a tricritical point (TCP) at small 
negative values of 7, but that was never confirmed in actual simulations at 
7 < 0. If there is indeed a 2 nd order transition, its properties have not been 
determined up today. 

On the other hand, both computational and data analysis techniques have 
progressed. This provides us with the possibility to perform a thorough FSS 
study of this model in a new context. Practically all other studies have dealt 
with the standard periodic boundary conditions, i.e. hypertorus topology for 
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the lattices. It has however been realized that there are nonlocal excitations 
in the system — closed monopole loops — that may extend over the whole 
lattice. Therefore the essentially local updating algorithms used for gauge 
theories, together with the boundary conditions, may affect thermalization 
properties. One expects, that in the thermodynamic limit the "continuum" 
properties of the system are independent of the global topology of the system, 
if this becomes locally flat. For these reasons it was suggested to simulate 
the model on lattices with spherelike topology [2], amounting to modified 
boundary conditions, such that closed loops are always homotopically equiv- 
alent to points. The spherelike topology allows the monopoles more freedom 
in their dynamics without changing the action. 

One generally expects that the thermodynamic properties of the bulk 
phase (defined by the behavior of the free energy per unit volume in the 
thermodynamic limit) are not affected by contributions which grow slower 
than the total volume. Boundary contributions are suppressed 0(1/ L) rel- 
ative to the leading term, curvature terms at least 0(1/ L 2 ), and therefore 
they should not change the critical exponents of the bulk phase. This does 
not necessarily hold for the ground state. E.g. at phase transitions of 1 st 
order the phase mixture may be different, depending on boundaries or even 
one individual spin, or due to an overall external field vanishing 0(1/ L D ). 
A similar influence may come from the curvature, even if thinned out over 
the volume in 0(1/L 2 ). However, if the manifold becomes locally flat in 
the thermodynamic limit the systems universal critical properties should be 
independent of the global topological structure. Otherwise we could hardly 
assume that we can do reliable physics on earth without definite knowledge 
about the topological details of the universe. 

Whereas in the original study the surface of a 5D cube was used, we 
here choose a discretization of the sphere, where the curvature is distributed 
more homogeneously over the lattice. Although the system is nonhomoge- 
neous on the scale of the lattice constant, it is homogeneous and isotropic 
on larger scales. As will be demonstrated here, in this 4D system with the 
topology of the surface of a 5D sphere we find no two-state signal on lattices 
with up to almost 20 4 points. Of course we cannot exclude the possibility, 
that a two-state signal reappears for even larger lattices. However, our FSS 
analysis leads to consistent results for a PT of 2 nd order for 7 < 0. 

The best measurement of the correlation length critical exponent v by 
means of the FSS behavior of the Fisher zero gives v = 0.365(8). Less 
precise FSS analyses of several other bulk observables are consistent with v 
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values in the interval v = 0.35 — 0.40. As we argue in t 3], due to rigorous 
dual relationships our results imply that also the Coulomb gas of monopole 
loops [4] and the noncompact U(l) Higgs model at large negative squared bare 
mass (frozen superconductor) [5] have a continuum limit described by the 
same non-Gaussian fixed point. Some further related models are discussed 
in 

The first order signal persists — also for the discussed spherelike geome- 
tries — at values 7 ~ 0.2. Since scaling may be garbled close to TCPs we 
concentrated on negative values of the double charge coupling 7 in our study. 
Scaling and FSS is expected to improve at larger distance from the 1 st order 
part of the PT line. Nevertheless, at 7 = the two-state signal is still absent 
and the scaling behavior is consistent with that found at 7 < 0. 

Let us add a remark on the extended action considered. Some time 
ago it was pointed out [7j, that, although the Wilson- and the heat kernel 
(Villain) action do have reflection positivity, some actions do not. Positivity 
is a sufficient, but not a necessary condition for unitarity; it guarantees the 
existence of a positive definite scalar product and the spectral condition, one 
of the formal conditions for the existence of a continuum limit field theory 
jH]. The actions with the parameter values 7 < considered here are not 
reflection positive. They share this property with other actions like e.g. the 
(Symanzik) improved actions. 

On the other hand, if reflection positivity holds on a part of a critical 
surface that is in the domain of attraction of a fixed point of some renormal- 
ization group transformation, we expect that it should be fulfilled everywhere 
in that domain. Unitarity violating states like ghosts should then decouple at 
large scale. We therefore find it justified to study the action near candidates 
for critical points even outside the region 7 > 0, where reflection positivity 
is respected on the scale of the lattice spacing. Unitarity at 7 < is also 
suggested by the observation that the regions with 7 < and 7 > are 
connected by the RG flows jUJ. 

Following a brief review of the situation in the U(l) pure gauge model we 
then introduce the spherical lattice in sect. El (Further technical details are 
given in the appendix.) In sect. Elwe present the Monte Carlo simulation and 
discuss the observables in some detail, including the expected FSS behavior. 
The results and their analysis are summarized in sect. EJ followed by our 
conclusions. 
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1.2 Situation of U(l) pure gauge studies 

U(l) is the most elementary Lie group that can be used to construct a quan- 
tum gauge field theory. Yet, when formulated on a lattice, the pure U(l) 
gauge theory turned out unexpectedly to be one of the most intriguing and 
less understood quantum gauge models. The awareness of this fact has de- 
veloped with an accumulation of the numerical experience. In this section 
we give a brief description of this historical development. We apologize for 
quoting only a subjectively chosen part of a much larger number of valuable 
papers. 

Since the introduction of lattice gauge theories by Wilson [TU] . the pure 
compact U(l) theory has been of interest as a theory with a rigorously estab- 
lished [TT] IT2"] PT between the confinement and the free charge (Coulomb) 
phases at zero temperature. One reason was the importance of topological 
excitations, the monopoles, for confinement, as manifested by their promi- 
nent role in this phase transition jU HBJ EH UE1 UH \H\ • Another purpose was 
to study it as a prototype example for applications of numerical methods of 
statistical physics in a lattice gauge theory, in particular an investigation of 
the continuum limit at the phase transition. However, the lesson has been 
that this phase transition provides no simple exercise. 

In the very first numerical investigations [U HHJ [THl EH1 EE], restricted 
to 7 = and small lattices, a behavior consistent with a 2 nd order phase 
transition at (3 ~ 1 was observed. But this order was questioned by the 
subsequent observation of a two-state signal on larger lattices [22 . Such a 
signal could imply that the phase transition at 7 = is actually of weak 1 st 
order, which would prevent taking a continuum limit there. The question 
was, and remained to be, whether this signal may be a finite-size effect. 

In the model with extended Wilson action ((TJ) it was found |2Sj that the 
confinement-Coulomb phase transition is clearly of 1 st order for 7 > 0.2, and 
weakens with decreasing 7. This suggested that the order of the transition 
changes when 7 is decreased, presumably at a TCP. 

The question at which value of 7 this happens turned out to be very 
difficult. First, even at large negative 7, a two-state signal has been observed, 
e.g. at 7 = —0.5 on the 8 4 lattice [21J. Second, TCPs are known to cause 
intricate finite-size effects |22]; easily mocking up a false order of the phase 
transition. 

In the hope to clarify the situation, an investigation of the strongly 1 st 
order part of the phase transition line at 7 > 0.2 was performed [21]. There 
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the latent heat Ae can be determined reliably even on moderately large 
lattices. Its independence on the lattice size was checked very carefully. The 
extrapolation of Ae to zero with decreasing 7 by means of the power law 

Aeoc( 7 -7 TC T u , (2) 

suggested that the order of the transition changes at the TCP with ^ TCP = 
—0.11(5), implying a 1 st order PT at 7 = 0. 

This extrapolation procedure is an attempt to control finite-size effects, 
but it uses the assumption that the power law behavior (J2J), which the data 
in the investigated region 7 = 0.2 — 0.5 are consistent with, indeed holds 
throughout the whole interval between 7 = 0.5 and 7 = r y TCP . This as- 
sumption has remained unverified. Another possible uncertainty in [21] was 
the determination of Ae at a strong 1 st order transition, without the more 
advanced methods of investigation of such transitions [26^ . 

Monte Carlo RG (MCRG) studies [2312313121 did not confirm this po- 
sition of the TCP at negative 7. Of course, also the MCRG approach suffers 
from ambiguities due to a small number of RG steps and a restricted number 
of couplings considered. Therefore also these studies remained inconclusive 
about the order of the transition around 7 = in the thermodynamic limit 
(although they all observed clear two-state signals). 

In spite of this, numerous attempts to determine the critical exponent 
v provided roughly consistent values in the range v ~ 0.28 — 0.42. These 
studies used various methods: the analytic calculations [33 [HI], the FSS 
analysis [13 E3 [23 E2| , the scaling of the string tension [231331 131 133, and 
the MCRG method [23 [23 13 133 123 • Three actions: Wilson-, extended 
Wilson-, and Villain-type have been used. 

This suggested that the pure compact U(l) lattice theory might have an 
interesting continuum limit at the confinement- Coulomb PT, pondered, e.g., 
in [3 EH] • However, the two-state signal, observed on finite lattices even for 
7 < [21] as well as for the Villain action [33) hindered the investigations 
of this possibility. Even if this signal is only a finite-size effect and the tran- 
sition in the infinite volume limit is genuinely of 2 nd order, it represents a 
serious impediment for a precise FSS analysis or MCRG studies. Because 
of this the investigation of the pure compact U(l) gauge theory lost its mo- 
mentum. Until now there is no established 2 nd order PT with an undisputed 
determination of critical indices in this model. 

All the above mentioned numerical work has been performed on 4D 
toroidal lattices. Recently, following earlier suggestions [33 E3 E] , two of 
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the present authors speculated that the two-state signal at 7 < may be 
related to monopole loops winding around the toroidal lattice, and trapped 
in simulations with local update algorithms [2 j. They used the 4D surface 
of a 5D cubic lattice instead of the torus, and observed that at 7 = the 
two-state signal vanishes on lattices of all investigated sizes. 

Choosing a spherelike topology provides a way to allow the monopoles 
more freedom in their dynamics without affecting them locally by changing 
the action. We consider this as preferable to adding terms to the action that 
forbid or restrict the occurrence of monopoles produces 0(L D ) contributions 
to the total free energy and thus changes the bulk properties of the system 
I4()j. In that case one explores the phase diagram in different regions of 
the space of couplings and the position of the phase transition in j3 moves to 
different values, depending on the extra couplings. None of these studies has 
led to phase transitions of 2 nd order, though. 

However, the cause for the two-state signal on the toroidal lattice, and 
of its vanishing on a lattice with the topology of a sphere, is not yet fully 
understood. Possibly the trivial 1 st homotopy group of such a lattice allows a 
smooth vanishing of winding monopole loops in simulations. But some other 
recent results do not seem to support this interpretation |4TH |4*H 14*2*]. 

On the other hand, for a study of the continuum limit on lattices with 
the topology of a sphere, a complete understanding of the dynamics of the 
two-state signal on a torus is not really necessary. What is required is a con- 
struction of a spherelike lattice which is homogeneous, in order to avoid the 
possibly related problems with the FSS analysis, encountered in j2]. Achiev- 
ing that in this paper, we hope to give a new momentum to the investigation 
of the confinement-Coulomb phase transition in the pure compact U(l) gauge 
theory on the lattice. A construction of the continuum limit appears now to 
be feasible. 

2 Spherelike lattices and lattice geometry 

In an attempt to formulate the theory without modification of the locally 
defined plaquette action and without forbidding or hindering the dynamic 
evolution of monopole loops, the lattice topology was modified. The usual 
periodic (or antiperiodic) boundary conditions correspond to the topology of 
a 4D torus T 4 with the first homotopy group Z 4 . Closed loops (or networks 
of loops) cannot necessarily be contracted to a point and the corresponding 
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ground state may be classified accordingly. The original motivation to divert 
from this lattice structure was to study the possible influence of this property 
on the dynamics of the phase transition. 

In [2] it was therefore suggested to simulate and study the model on 
a lattice with spherelike topology, in particular on the surface SH[N] of a 
5D hypercubic lattice of size N 5 . This lattice may be best visualized in 
analogy to the 2D surface of a 3D cubic lattice. It may also be considered 
as a collection of 10 hypercubic lattices of size iV 4 , glued together at their 
boundaries. This implies that one can expect the same critical coupling in 
the thermodynamic limit as for the usual torus. This was indeed verified in 
the Monte Carlo calculations |2j- Details and parameters of the geometry 
are listed in the appendix. 

The so-defined lattice is locally flat, except at certain plaquettes ((D- 
2)-dimensional elements), where the curvature is concentrated, a well-known 
property of Regge skeletons. The unusual features include plaquettes border- 
ing only three 3D cubes (instead of the usual 4), links bordering less than 6 
plaquettes and sites with less than 8 links. These curvaturelike contributions 
as we might call them in the absence of a strict theory in 4D are suppressed 
0(1/N 2 ) relative to the leading terms in the action. 

In an attempt to distribute these local inhomogeneities more uniformly 
over the lattice we introduced the "almost smooth" spherical lattice S[N]. 
In the construction we project sites, links and plaquettes of SH[N] (or its 
dual SH'[N]) onto the surface of a concentric 4D sphere and introduce weight 
factors similar to those used by [13] hi their study of random triangulated 
lattices, 

S = -J2 w P [(3 cos(e P ) +7cos(26 P )] 
p 

A' 

with w p = ——. (3) 
Ap 

Here, A P and A' p denote the areas of the corresponding plaquette and its 
dual, respectively, of the projected lattice. 

As discussed in j3H] in the situation of triangulated random lattices, one 
has to distribute the total integration volume over all contributions to the 
action, i.e. the plaquette terms in our case. This may be done with the help 
of the dual lattice, where to each site, link or plaquette there is an associated 
dual 4D cube, 3D cube or plaquette. The dual lattice sites in our situation are 
constructed from the barycenters of the 4D cubes that have been projected 
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to the sphere S 4 , followed by an adjusting projection of these points to the 
sphere. Further reasoning according to [13] leads to the weight factors wp in 
Q. This choice is not unique, but reproduces the usual continuum action 
F? v in the naive continuum limit g — > (j3 —>■ oo) and is thus homogeneous in 
this limit. We study the system at finite /?; there slight distortions from the 
regular spherical surface are possible. The value of (3 cr i t might be modified 
due to the weight factors and thus does not necessarily agree with that of 
the torus or SH. 

Technically we have to introduce some approximations. Usually the pla- 
quettes - constructed via the projection of the sites to the sphere - will not 
be flat. Plaquette areas are therefore determined from the sum of two tri- 
angles. Also for this reason the curvature in this formulation will not be 
perfectly uniformly distributed. In order to achieve this, we would have to 
rely on a triangulated lattice. This would imply a significant change of the 
action, which we wanted to avoid. On the other hand we expect these ef- 
fects to become irrelevant in a situation with a large correlation length. The 
consistency of the found FSS behavior justifies these simplifications. 

The connectivity properties of SH and S are identical. In the computer 
programs the geometry is implemented with tables and the weight factors 
wp are precalculated. During the development of the program and in the 
early stages of the analysis we also determined the monopole positions (on 
the dual lattice) and reproduced them graphically. We observed the expected 
properties: 

• The monopole loops were always closed; 

• the smallest loops had length 3 (corner plaquettes on the dual lattice); 

• they fluctuated freely, appearing and disappearing without noticeable 
correlation with positions close to corners. 

This also served as a check of the consistency of the connectivity tables. 
In our discussions we will refer to the effective lattice volume 



as the typical size quantity. For SH this would be just the number of sites, 
for S it is very close to this value. A length scale may be defined as 



u p 



(4) 



L 



V*. 



(5) 
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Table 1: Effective volumes for the studied lattices SH[N]. We also give the 
value of L — which would give the base length for a hypertorus lattice 
with the same volume. 



N 

1 V 




T, 


4 


825.1 


5.4 


5 


2576.6 


7.1 


6 


6268.1 


8.9 


7 


12986.9 


10.7 


8 


24064.1 


12.5 


9 


41074.6 


14.2 


10 


65837.3 


16.0 


12 


147113.8 


19.6 



Table [T] summarizes the effective volumes for the lattice sizes used in our 
study. 

For strictly asymptotic dependencies as they come up in FSS studies it 
is irrelevant, whether one uses iV or L. However, for moderately sized finite 
systems a suitable choice improves the approach to the asymptotic behavior. 
Let us mention in this context that in |H] different lattice geometries were 
compared and it was demonstrated, that the scaling curves show best agree- 
ment with each other, if one indeed uses L - the size derived from the total 
volume as opposed to the base length N - as size variable. In the present 
work we cannot compare with torus results, since for those the two-state 
signal obscures the measured values of the cumulants. 

3 Simulation methods and statistics 

In an earlier study [2] we found that for lattice type SH there are two-state 
signals at the pseudocritical (3 for 7 = 0.2, but no such indications at 7 = 0. 
For this reason we now studied the action (J3J) of the spherical lattice S at 
the values 7 = 0., —0.2, —0.5. Preliminary results have been presented in 

gaum. 
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3.1 Updating and measuring 

We have worked with lattices S[N] for N ranging between 4 and 12. The 
couplings were chosen in the immediate neighborhood of the pseudocritical 
values of (3. For the analysis we determined the histograms of the weighted 
sum of plaquette values 

E = wp cos Op. (6) 
p 

Note that this is not the total energy, but just the part corresponding to the 
coupling parameter /3; it coincides with the total plaquette energy for the 
Wilson action, 7 = 0. Any scaling- or two-state signal should be observable 
in that quantity. We also define the density 

e = E/^2w P = E/(6V). (7) 
p 

We combined the various histograms (for fixed 7 but different /3) with help 
of the Ferrenberg-Swendsen multihistogram reweighting technique |47j . For 
each 7 we thereby construct the density of states p(E; 7). The representation 

Z(/3, 7 ) = 5>(£; 7 )exp(-/3£) (8) 

E 

allows us to determine various observables for continuous values of f3 through 
(E n ) = j^- ) J2 P (E; 1 )E n exp(-l3E). (9) 

Since we never observed two-state signals we did not implement multicanon- 
ical updating j2H|- We used a 3-hit Metropolis update; for 7 = we included 
an additional overrelaxation step (the autocorrelation length decreased by 
a factor of about 2). For each lattice size at each 7 we typically accumu- 
lated O(10 6 ) updates, which is between 10 3 and 10 4 times the integrated 
autocorrelation length for the observable E (cf. Table |2J). 

The histograms had up to 10 4 bins in order to exclude any possible influ- 
ence due to binning. In fact, by rebinning we found no changes within single 
precision down to 0(500) bins. Due to the fine binning the raw histograms 
have a noisy appearance, which is irrelevant for the Ferrenberg-Swendsen 
analysis and for the representation (jHJ). For the plots we use rebinned ver- 
sions with only 500 bins maximum. 

The total CPU-time spent for the calculations on workstations and on a 
Cray-YMP sums up to 6800 hours in Cray-YMP units. 
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Table 2: Statistics of the data for the studied values of 7: Lattice base length 
N, total number n of configurations in multiples of 10 6 , range of /3-values, 
maximal r int values (from a fit to all beta values, as discussed in the text; 
errors are typically 10% of the values). 
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A T 


/1 n6 
72/10 


/?-range 


T~int,max 


0. 


4 


1.1 


0.98- 


-1.025 


12 




6 


1.0 


1- 


-1.025 


56 




8 


1.1 


1.0125- 


-1.0275 


150 




10 


1.47 


1.01- 


-1.025 


304 


-0.2 


4 


1.1 


1.07- 


-1.2 


24 




5 


1.8 


1.15- 


-1.175 


63 




6 


1.6 


1.13- 


-1.21 


116 




7 


1.1 


1.14- 


-1.175 


188 




8 


1.6 


1.155- 


-1.175 


272 




9 


1.1 


1.1655- 


-1.1715 


367 




10 


1.0 


1.1635- 


-1.1715 


583 


-0.5 


4 


4.8 


1.3- 


-1.65 


37 




5 


1.8 


1.38- 


-1.47 


76 




6 


2.1 


1.35- 


-1.5 


149 




7 


1.9 


1.402- 


-1.452 


311 




8 


1.5 


1.4- 


-1.455 


332 




9 


1.8 


1.42- 


-1.442 


480 




10 


1.55 


1.42- 


-1.442 


473 




12 


1.6 


1.43- 


-1.44 


1565 



3.2 Observables and FSS 
3.2.1 Cumulants 

We determined the 2 nd and 4 th order cumulants of the observable E. Due 
to the analogy to the internal energy (identical to E only for 7 = 0) we call 
the second order cumulant the specific- heat. The specific-heat, the Challa- 
Landau-Binder cumulant [IB] and another 4^ order cumulant suggested by 
Binder (cf. [49] and jSOj) are defined through 

cv(P,L) = ^((E-iE)) 2 ), (10) 
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v (3D - 1 ~ ^ 2 » 2 > fin 

U(3L) - « E ~^ (12) 

The positions and values of their respective extrema are used for the FSS 
analysis. 

From the usual scaling hypothesis [HH E21 EBJ E~Q EH] one expects for the 
singular part of the free energy density the scaling behavior 

f(r,L) = L-^ D f(rL^,l), (13) 

where r = (1 — 3/3 c ) denotes the reduced coupling and L is a length scale 
parameter. From this one derives the scaling behavior of the cumulants. At 
a 2 nd order phase transition we expect (for D=4 and a > 0) 

C V ,ma X (L) ^ L a '\ (14) 
V CL B,min(L) ~ L^" 4 , (15) 

Ut, min (L) ~ 0(1)+0(L-^), (16) 
P e (L)-/3 c ~ L- A . (17) 

For a = there are logarithmic terms. The asymptotic value of U^ m i n 
depends on the details of the distribution density p(E) and is 3 for a Gaussian 
distribution. Mean field values are v — 1/2 and with Josephson's law a = 
2 - Dv = 0. 

We denote by /3 C (L) definitions for pseudocritical points like the positions 
of the extrema in the cumulants. The so-called shift-exponent A is for many 
models equal to 1/V, but not necessarily so in general; such an identity is 
not a necessary result of FSS (cf. the discussion in We return to this 

issue later. Furthermore, a priori we know nothing about the absolute size 
of the multiplicative coefficients in the scaling formulas. They depend on the 
details of the action, the lattice geometry and the topology [53] . 

For 1 st order transitions one expects the FSS behavior 

L~ D c Vtmax -> \{e -e d ) 2 , (18) 
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V C LB,min - ^ + 0(L~ D ), (19) 

U 4 , min -> 1 + 0(L~ D ), (20) 
3 C (L)-3 C = 0(L- D ). (21) 
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Here e Q and denote the discontinuous values of the energy density at the 
PT point. 

As discussed, in the considered lattice geometry there are inhomogeneities 
in the sense, that the coordination numbers of some sites, links and plaquettes 
deviate from the usual torus numbers. For the SH lattices these may be 
considered as lattice inhomogeneities. Their contribution to the total free 
energy is suppressed 0(1/N 2 ) = 0(V~ 2 / D ). In our smoothed version of 
that lattice: S, the inhomogeneous contribution should be smaller. There is 
however still a possible contribution of the total curvature to the free energy, 
which is suppressed with the same order (cf. also the discussion for the 
2D Ising model 01]). Thus — in principle - we also may expect "surface" 
corrections of 0(V~ 2 / D ) in all FSS relations. Indeed such contributions have 
been observed for the SH lattices [2|. It turns out that they are much smaller 
in our present study, in fact too small to study them. 

3.2.2 Fisher zeros 

Equation (JHJ) defines implicitly an analytic continuation to complex values of 
(3 not too far away from the real axis. Therefore it is possible to determine 
the nearby zeros of the partition function [HE] in the complex /3-plane, the 
so-called Fisher zeros j^Zj (cf. [HE]). 

One should add a warning concerning technical aspects. The histograms 
are binned, having both, upper and lower limits for E max and E min as well as 
a bin size A = (E max — E min ) x 10~ 4 . The representation (JHJ) for /3 = (3r + %(3i 
therefore is a discrete Fouriertransformation. It will induce a periodicity in f3i 
due to the bin size and an effective grid with grid spacing 27i/(E max — E min ) 
(although the values of the partition function Z are well-defined even between 
the grid points, they carry no additional information). 

Usually the distribution is similar to a Gaussian distribution; let us for 
the sake of the argument assume such a form 

p(E) exp (-P R E) ~ exp (-c{E - E ) 2 ) . (22) 

From (JHJ) one then expects an oscillatory behavior of Z proportional to 
exp (if3iE ). This is indeed observed in the calculation. In the search for 
partition function zeros one starts with an identification of sign changes of 
Im Z and Re Z. The rapidly oscillating phase factor may confuse the pattern 
and one has to work with a very fine resolution and to carefully combine 
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the sign-change analysis with a search in \Z\. Also the grid structure may 
interfere with these oscillations and one has to proceed with care. 

So real and imaginary parts of the closest Fisher zeros provide further 
(even) observables. In particular the imaginary part of the zero closest to 
the real axis provides a high quality estimator for the critical exponent v 
(cf. [HE] for a recent high statistics study of the 4D Ising model, where it 
was possible to identify the logarithmic corrections to scaling on basis of the 
Lee- Yang zeros |56j). As will be demonstrated below this quantity appears to 
have small corrections to the leading FSS behavior in our environment; this 
observation is analogous to other recent investigations of spherelike lattices 

in 2D isamj. 

From the scaling arguments for the free energy we expect 

|2o(X) — P c \ — L" 1 ^ . (23) 

This provides an upper bound for the real and imaginary part of Zq, in 
particular 

lmz (L) ~ 0(L- 1/u ), (24) 

Rez (L)-f3 c ~ 0{L- l ' v ). (25) 

Although in some cases the angle, under which the zeros approach the real 
axis (defined as the angle of a line connecting the two closest zeros) is known 
(e.g. 7r/2 for the 2D Ising model in the Onsager solution [2Zj, 7r/4 in the 
mean field solution for the 4D 4 -model [BT] , both on cubic lattices with torus 
topology) there is no FSS theory for this angle of approach with regard to the 
size L. Depending on details of the model, the geometry and the topology, 
the real part — which by analogy to the cumulants we call a pseudocritical 
value — may approach the asymptotic value faster, i.e. with a shift exponent 
A larger than 1/v [53] . Such a behavior has been observed in a recent study 
of the 2D Ising model jH]. 

We also mention here, that the position of the closest Fisher zero is re- 
lated to the peak position and value of the specific-heat. Since the partition 
function may be expressed by the Vieta-product of all its zeros {zi}, the 
specific-heat is proportional to 

Vwh? (26) 
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and therefore the leading contribution to V x cy near the PT is propor- 
tional to (Im^o)" 2 . The peak position is in leading order given by Rezo- Of 
course there are further contributions due to the other zeros and a possible 
background from an entire function. 

Also these observables may in principle exhibit corrections to FSS due to 
curvature and topology as discussed above. 



4 Results and data analysis 

4.1 Autocorrelation and error analysis 

For all individual runs we determined the integrated autocorrelation for the 
observable E, 

(Here the index indicates the n-th configuration measured in the Markov pro- 
cess.) The inverse value provided us with a weight factor of the corresponding 
data sample in the multihistogram analysis. 

For the maximum values of the autocorrelation lengths we produced two 
sets of values. One was the maximal observed values of T int ^ for all samples 
for the given lattice size. The other resulted from a fit to the values of 
Tint{P>) t° a peak shaped curve which has its peak position where the specific- 
heat (see the discussion below of the analysis for the cumulants) assumes its 
maximum. In the subsequent analysis we discuss only the results due to the 
first set. The second set led to similar results. 

The maximum values of Tint s (cf- Table |2J) increase from values of ~ 12 
for 7 = 0, L = 4 up to ~ 1600 for 7 = —0.5, L = 12. This demonstrates the 
necessity to work with large samples of several 10 5 configurations for each 
value of /3, at least for the large lattices. 

At second order PTs the maximum values of the integrated autocorre- 
lation time provides an estimate for the corresponding dynamical critical 
exponent zs through 

r inM; ~mm(L,0 ZE (28) 

(where £ denotes the correlation length). At 1 st order transitions one expects 
that the autocorrelation length grows exponentially ~ exp (cL D_1 ). In FigQ 
a log-log plot shows that the size dependence is indeed compatible with 
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Figure 1: The maximum values of the integrated autocorrelation length for 
the energy observable together with the fit to an exponential dependence on 
the lattice size. 

We assume, that the peak values correspond to the point where 

TintAPpeak) = c(*y)L' E . (29) 

A simultaneous fit to all three data sets (for the three values of 7) gives for 
the dynamical critical exponent ze = 2.85(6). The coefficients grow from 
c(0) = 0.12 up to c(— 0.5) = 0.32. The results for 7 = were obtained with 
an additional overrelaxation step in the Metropolis updating. Although the 
absolute value of the autocorrelation lengths decreased by a factor of about 
2, the dynamical critical exponent appears not to be affected. 

That value is substantially larger than the value z = 2 expected for the 
random walk dynamics of local algorithms and demonstrated for Gaussian 
models. This behavior is indicative of a more complex dynamics than it is 
usually anticipated for systems with local excitations. The non-locality of 
the monopole loops may be responsible for the observed effect. On the other 
hand, we may not yet be asymptotic and the determination of reliable values 
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for this exponent is notoriously difficult. 

As a consistency check we also determined autocorrelation times from a 
fit to an exponential decay and from blocking analysis. The resulting values 
were typically proportional to those discussed above, although less reliable, 
i.e. with large statistical fluctuations. The exponential autocorrelation time 
and its dynamical critical exponent are upper bounds to the integrated au- 
tocorrelation time (cf. [62 J. 

The statistical errors for all our raw data (i.e. positions and values of 
cumulant extrema and positions of the Fisher zeros) were determined with 
the jackknife method. From the original set of values E for each configuration 
we chose 10 different subsets by omitting 10% of the numbers, providing 
10 histograms. The Ferrenberg-Swendsen analysis then was repeated for 
all these subhistograms and parameters for the cumulants (peak positions, 
values, Fisher zeros) were determined. The distribution of these numbers 
defined the errors according to the jackknife procedure. The central values 
were taken from the analysis of the complete data. The fits were performed 
using these central values and errors. 

The simulations on the CRAY-YMP have been performed employing a 
vectorized version of the shift-register random number generator, which in 
its actual implementation uses XOR operations in between the i and i + 103 
element to generate the i + 205 element of the sequence. For the programs 
on the workstations we used a corrected version of RCARRY [OS] based on 
the "subtract-and-borrow" version of a lagged Fibonacci algorithm. 

4.2 Results: data and fits 

We analyzed the final numbers for the pseudocritical points (the extrema 
positions of the cumulants and the real part of the position of the closest 
Fisher zero), the extrema values of the cumulants and the imaginary part of 
the Fisher zero. The fits were performed both, for all lattices sizes and for a 
subset of lattices with N > 6, in order to estimate to which amount we see 
asymptotic behavior. 

4.2.1 Histograms 

From the combination of histograms determined for different values of (3 
according to the Ferrenberg-Swendsen technique we obtain the distribution 
densities /?(.£; 7) in (J5|). A necessary condition for the effectiveness of the 
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Figure 2: Pseudo-histograms as function of e at 7 = —0.2 and for lattices 
size N =4, 5, 6, 7, 8, 9, 10 (from left to right) at the respective peak positions 
of the specific-heat (in Table El) and normalized to unit maximum value. 

approach is sufficient overlap between the individual histograms. From the 
densities one may construct pseudo-histograms (or reweighted histograms) 
at arbitrary values of (3 (which should be in the domain covered by the 
individual histograms) , 

7) = />(£; 7) exp (-/?£). (30) 

These interpolate the individual histograms but they also bring together and 
represent all histogram data. 

In Fig. El we plot the pseudo-histograms for 7 = —0.2 for all lattice sizes 
studied and determined at the peak positions of the specific- heat. They are 
normalized to unity at their respective maxima. No double peak structure 
is observable. In Fig. El the pseudo-histograms for the largest lattice sizes 
are plotted for all three values of 7. (An individual histogram at 7 = is 
shown in jHj.) Again, there is no indication of a discontinuity signal. Such 
an observation was made already in the study of the SH lattice (at 7 = 0) in 
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Figure 3: Pseudo-histograms vs e for lattices (from left to right) S [12] 
( 7 = -0.5), S [10] (7 = -0.2) and S [10] (7 = 0) at their respective peak 
positions of the specific-heat and normalized to unit maximum. 

® 

Also the individual histograms show no two-peak structure. Actually, 
if the statistics is small, spurious signals may appear, but with increasing 
statistics they always vanished. (They also were not at consistent positions.) 
Finally we mention, that there are studies at established but weak first order 
transitions (the 2D 5-state Potts model) on similar spherelike lattices, where 
a two-state signal has been observed clearly [46] . Therefore we find no in- 
dication, that the particular kind of lattice studied here has a tendency to 
smear out weak two-state signals. 

A two-peak distribution is an indicator of a possible 1 st order transition. 
However, in order to establish this order one should also find further signals 
for coexistence of phases like FSS consistent with v = 1/D = 1/4, correct 
scaling of the minimum between the peaks (suppressed due to the 3D in- 
terface) in the distribution and tunneling probability oc exp(— 2crL 3 ). Up 
to now no consistent observations of that kind have been made in the U(l) 
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theory for the toroidal lattices, where one finds two-state signals. 

Within the scope of lattice sizes studied here we are therefore led to 
assume a 2 nd order behavior for 7 < 0. The subsequent scaling analysis leads 
to results fully supporting this assumption. 

In Tables H3 and 0] we summarize our results for the extrema values and 
positions of the cumulants and of the positions of the closest Fisher zeros. 
The analysis of these data is discussed in the subsequent sections. 

Table 3: Extrema positions of the cumulants and the real part of the positions 
of the c losest Fisher zeros. 



7 


N 


P(pv) 


P{V C lb) 


(3{U A ) 


Re(^o) 


n 


4 


1 0027(3") 


9990(4) 


1 0051 (3) 


1 (1047(4^ 

X . uu^ 1 It:/ 




6 


1 0151(2) 


1 0148(1) 


1 0156(2) 


1 0156(2) 




8 


1.0179(1) 


1 0179(1) 


1 0182(1) 

X . W X I X J 


1 0182(1) 




10 


1.0183(1) 


1.0183(1) 


1.0185(1) 


1.0185(1) 


-0.2 


4 


1.1473(6) 


1.1422(5) 


1.1512(11) 


1.1514(12) 




5 


1.1588(4) 


1.1574(5) 


1.1607(5) 


1.1608(5) 




6 


1.1640(3) 


1.1634(3) 


1.1652(7) 


1.1650(12) 




7 


1.1664(4) 


1.1662(4) 


1.1675(4) 


1.1677(2) 




8 


1.1681(1) 


1.1680(1) 


1.1685(3) 


1.1684(3) 




9 


1.1688(1) 


1.1687(1) 


1.1690(1) 


1.1690(1) 




10 


1.1695(1) 


1.1695(1) 


1.1698(2) 


1.1697(2) 


-0.5 


4 


1.4067(7) 


1.3987(10) 


1.4126(15) 


1.4070(42) 




5 


1.4202(7) 


1.4177(8) 


1.4239(13) 


1.4246(17) 




6 


1.4270(6) 


1.4262(6) 


1.4291(5) 


1.4292(7) 




7 


1.4307(4) 


1.4304(4) 


1.4320(6) 


1.4318(5) 




8 


1.4325(2) 


1.4324(2) 


1.4328(3) 


1.4328(3) 




9 


1.4340(7) 


1.4339(8) 


1.4354(22) 


1.4353(4) 




10 


1.4346(2) 


1.4345(2) 


1.4349(2) 


1.4349(2) 




12 


1.4359(6) 


1.4359(6) 


1.4366(1) 


1.4365(1) 



4.2.2 Fisher zeros 

The results for the imaginary part of the position z of the Fisher zero closest 
to the (3 axis are given in Table El Although we tried fits including further 
background contributions it turned out that the form (|2^j) is sufficient. 
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Table 4: Extrema values of the cumulants and the imaginary part of the 
positions of the closest Fisher zeros. 



7 


N 


cv 


VqlB 


u 4 


Im(z ) 





4 


1.85(1) 


-0.142(l)E-02 


2.77(1) 


0.0300(5) 




6 


3.93(5) 


-0.362(5)E-03 


2.67(2) 


0.0066(2) 




8 


6.75(14) 


-0.157(3)E-03 


2.61(3) 


0.0024(1) 




10 


9.47(23) 


-0.792(20)E-04 


2.65(2) 


0.0013(1) 


-0.2 


4 


1.22(1) 


-0.982(6)E-03 


2.83(1) 


0.0388(10) 




5 


1.63(1) 


-0.395(3)E-03 


2.81(1) 


0.0185(4) 




6 


2.07(3) 


-0.201(3)E-03 


2.82(2) 


0.0112(14) 




7 


2.54(5) 


-0.117(2)E-03 


2.80(4) 


0.0063(3) 




8 


3.08(4) 


-0.755(9)E-04 


2.80(2) 


0.0043(2) 




9 


3.54(9) 


-0.503(12)E-04 


2.75(5) 


0.0029(2) 




10 


4.22(13) 


-0.371(ll)E-04 


2.71(8) 


0.0020(2) 


-0.5 


4 


0.76(1) 


-0.647(3)E-03 


2.89(1) 


0.0578(32) 




5 


0.95(1) 


-0.246(3)E-03 


2.88(1) 


0.0271(11) 




6 


1.16(1) 


-0.121(l)E-03 


2.85(2) 


0.0150(11) 




7 


1.36(1) 


-0.673(8)E-04 


2.83(1) 


0.0087(2) 




8 


1.62(6) 


-0.427(16)E-04 


2.72(6) 


0.0053(4) 




9 


1.69(5) 


-0.258(8)E-04 


2.87(6) 


0.0046(5) 




10 


1.98(8) 


-0.188(8)E-04 


2.78(5) 


0.0030(2) 




12 


2.26(12) 


-0.946(49)E-05 


2.78(12) 


0.0019(2) 



Table 5: Results for v from individual fits to Imzo according to (J2IJ). 



7 


V 


X 2 /d.f. 





0.345(3) 


4.7 


-0.2 


0.378(7) 


0.3 


-0.5 


0.368(8) 


0.8 
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Figure 4: Log-log Plot for Imzo vs L with individual fits for each 7. At 
7 = the full line denotes our fit to all N, the dotted line a fit to data for 
N > 6. For 7 < no visible difference would be seen. 

In Fig. |3] we show individual fits for each 7. For 7 = we distinguish two 
lines: the fit to all N and one to data for N > 6. In particular for 7 = 
the N = 4 data seems to be outside the overall behavior, indicating that at 
this lattice size the asymptotic behavior is not yet seen. According to our 
interpretation, we expect the value 7 = to be closest to a tricritical point, 
which may explain the larger deviations as compared to the other values of 

7- 

A joint fit to the data for all three 7-values with universal v but indi- 
vidual proportionality factors gives v = 0.354(3) at a x 2 /d.f. value of 2.7; 
including only data with Af > 5 we obtain v = 0.368(5) (x^/d.f. = 0.98). 
Finally if we restrict the fit to the data with iV > 6 we find v = 0.365(8) at 
a {x 2 /d.f =1.05). This last fit we consider to be the most reliable determi- 
nation of v (a corresponding plot may be found in [3]). 

It is interesting to compare the absolute positions of the zeros for different 
values of 7 in Fig. El We find that the zeros are generally closer to the real 
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axis for 7 closer to 0. This indicates, that asymptotic scaling sets in somewhat 
later (at larger lattices) at more negative values of 7. This correlates with the 
peak values of the specific-heat, as will be discussed below in the discussion 
of the cumulants. 

The results for the real parts of the Fisher zero positions will be discussed 
together with the pseudocritical values. 
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Figure 5: Positions of the closest Fisher zeros for all lattice sizes and all 7. 



4.2.3 Cumulant values 

Figures El and [3 exhibit (e) and cy in the pseudocritical range, and Fig. |H1 
gives an example for the behavior of Vclb- The inserts in Fig. [^demonstrate, 
that the peak values of the specific-heat grow slower than the volume and 
that cy /V approaches zero in the thermodynamic limit, indicating a 2 nd 
order PT. 
Our ansatz 

c v ,ma X (L) = a(j) + 6( 7 )L 2 /^W- 4 (31) 
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Figure 6: Values of (e) vs (3 around the pseudocritical points for the three 
values of 7 and for all lattice sizes: (a) 7 = 0, N = 4, 6, 8, 10; (b) 7 = —0.2, 
iV = 4, 5, 6, 7, 8, 9, 10; (c) 7 = -0.5, iV = 4, 5, 6, 7, 8, 9, 10, 12. 

for the scaling behavior (|T4*|) (with Josephson's law relating a with v) al- 
lows for a background constant. We performed various fits restricting the 
coefficients in different ways. It turned out, that one should not omit the 
background constant 0(7). If one does, then the fits become size dependent 
and have worse x 2 if a U lattice sizes are included and better if one omits 
the small lattices. We therefore allow for such a background parameter and 
include all lattices sizes in the fits. 

If we leave a, b, and v 7— dependent we get consistent results with v 
varying between 0.361 and 0.404 (cf. TableEJ)- If we enforce a 7-independent 
value of v we find v = 0.378(4) and the reasonable \ 2 = 2.2 but different 
values for the background parameters. 
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Figure 7: Values of cy vs (3 around the pseudocritical points for the considered 
values of 7 and for all lattice sizes: (a) 7 = 0, N = 4, 6, 8, 10; (b) 7 = —0.2, 
N = 4,5,6,7,8,9,10; (c) 7 = -0.5, N = 4,5,6,7,8,9,10,12. The inserts 
exhibit the peak values cy /V vs 1/V demonstrating their approach towards 
for V — » 00. 



Table 6: Results of the fit to cy according to (pH|). 



7 


V 


0(7) 


6(7) 


X 2 /d.f. 





0.361(6) 


0.07(18) 


0.136(34) 


2.3 


-0.2 


0.374(6) 


0.35(9) 


0.090(20) 


0.3 


-0.5 


0.404(9) 


0.10(10) 


0.132(44) 


1.0 
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Figure 8: Vclb vs /3 for all lattices sizes studied at 7 = —0.2. For 1 st order 
PTs the values at the minima should asymptotically approach a nonzero 
constant. 

Fig. Elis a log- log plot for cy and the results of the fits (Table EJ). The 
increase of the value for v with decreasing 7 indicates that the behavior of 
cy is not yet asymptotic. We observed already in the discussion of the Fisher 
zero that scaling appears to be retarded towards more negative values of 7. 
Below (in I4.2.5j) we try to correct for this fact by introducing a phenomeno- 
logical scaling variable. Indeed we find a consistent scaling behavior of the 
specific-heat maximum corresponding to a value of v as determined in sect. 
I^2~2l 

For the CLB-cumulant we found that a fit to the FSS behavior in the 
form 

V C LB,min(L) = (o( 7 ) + b^L 2 '^) L~ 4 (32) 

- in the spirit of the correction term in the specific-heat (|3*Tjl - appears to 
be suitable. Fig. [TO and Tabled show our result. The consistency with the 
results for cy is remarkable. 

The values of Vclb clearly tend to vanish in the thermodynamic limit 
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Figure 9: This is a log-log plot of specific-heat cy vs L together with the fit 
results to form (j31|) and parameter values from Table El 

as expected for 2 nd order transitions. It was unexpected, that the scal- 
ing analysis led to sensible results in good agreement with the results for 
the specific- heat. (Note that the CLB-cumulant is a 4 th order moment and 
therefore in principle much more error-prune. For the - in comparison to 
spin model simulations - low statistics one cannot put too much confidence 
in this quantity.) 



Table 7: Results for the fit of VcLB,min according (j3*2Jl . 



7 


V 


0(7) 


6( 7 ) 


X 2 /d.f. 





0.361(6) 


-0.23(1) 


-0.071(19) 


2.3 


-0.2 


0.365(6) 


-0.41(5) 


-0.034(9) 


0.3 


-0.5 


0.396(9) 


-0.22(5) 


-0.054(19) 


1.1 



If we omit lattice sizes N < 6 the \ 2 improves, but the values of v do not 
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Figure 10: A plot of In (—Vclb) vs InL for all three value of 7, comparing 
with the fit according to (J32*j) . 

change much. Like for the specific-heat we notice an increase of v fit-values 
with more negative 7 which we interpret as due to the retarded FSS. 

The data for U4 show too little size dependence (or have too large errors) 
to produce a trustworthy fit to the expected leading scaling behavior (fTflj). 

U 4imm (L) = a( 1 ) + b( 1 )L 4 - 2 ^. (33) 

A joint fit to all data (x 2 /d.f. =1.1) gives v = 0.35(3) and values of a = 
2.60(3), 2.79(2), 2.81(2) which are, however, clearly different from the value 
1 expected at a 1 st order PT. 

4.2.4 Pseudocritical values 

Let us denote our four definitions for pseudocritical values by (3^ (L) (where 
% = 1...4 stands for the peak positions of cy, Vclb and U4 and Rezo, 
respectively). In the fits we allow for the form 

(5f{L) =(3 c + aiL- x . (34) 
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For each 7 we fit simultaneously to all types % for a unique j3 c and A but 
individual values a^. We find that allowing for another term 0(L~ 2 ) — as 
it is motivated from the possible contribution of the curvature or lattice 
inhomogeneities and as it seemed to be necessary for the analysis of the SH 
results in |2] — does not improve the \ 2 significantly 

Table 8: Results of the fits to the finite-size dependence of the four definitions 
of pseudo critical points according to (J53)l . 



7 


Pc 


1/A 


X 2 /d.f. 


0. 


1.0190(1) 


0.321(7) 


5.8 


-0.2 


1.1709(2) 


0.386(10) 


1.6 


-0.5 


1.4381(1) 


0.472(12) 


2.5 



TableElgives the fit values for the pseudocritical points and 1/A. It is not 
generally true, that A = \/v (cf. the discussion in [33]) and indeed a recent 
study indicated a different value for the 2D Ising model on spherelike lattices 
|44j . Accepting this caveat we still find numbers of similar size. If we allow 
for a correction term 0(L~ 2 ) due to the background curvature of our lattices 
and fix A = 1/0.37 (i.e. at a value 1/u suggested from the other data) the fit 
is of comparable quality with compatible values for P c (l) and the fit curves 
in the plots are indistinguishable by eye. 

Altogether the errors on the pseudocritical points are larger but the fits 
are not very satisfying (cf. Fig. ^2 as an example; the data and fits for the 
other 7- values look similar). The value of A is not stringently determined by 
the data (or the theory). 

4.2.5 Scaling consistency 

At finite lattices there are always corrections to FSS, depending on size, 
geometry, topology and of course details of the action and the observables. 
Since lmz gives the cleanest FSS signal, we use it as a phenomenological 
scaling variable 

x = lmz (35) 

and study in this section the other observables as functions of x. (Notice that 
x is defined from the data!) That provides us on one hand with a consistency 
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Figure 11: Fits to the data (for 7 = —0.2) for the pseudocritical points 
(according to the four different definitions described). We use the abscissa 
variable L~ l l u (for the preferred value v = 0.365) in order to emphasize the 
nontrivial dependence. 

check for our results. On the other hand this assumption allows us to bring 
together and combine results from different values of 7. 

Assuming the FSS relation x ~ Lr x l v (i.e. for Imzo) we expect e.g. for 
the specific-heat the behavior (cf. 1)140) 

c v , max (L) + L 2 '^ ~ a + fa^L" 4 , (36) 

where we - as discussed - allow for the additive constant to represent the 
unknown background. Fig. El shows the peak values of the specific-heat for 
the data vs x~ 2 L~ 4 together with the result of a linear fit to the data for N > 
5. For this fit we had to assume different values of the background constant 
a = 0.22,0.55,0.40 for different 7 = 0,-0.2,-0.5. The plot demonstrates 
the consistency of the FSS of the specific-heat values with that of Imzo. This 
is not surprising, since these quantities have a close relationship as discussed 
in sect. El 
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Figure 12: A plot of the peak values of the specific-heat (all 7, all lattices 
sizes) vs l/(x 2 L 4 ) (x is the phenomenological scaling variable defined in ([35)1: 
the corresponding error bars are shown as well. The line represents a linear 
fit. 

4.2.6 Summary of the fit results 

The cleanest and most consistent results on FSS come from the imaginary 
part of the closest Fisher zero positions and suggest a value v = 0.365(8). 

The peak values of the various cumulants are generally consistent with 
v ~ 0.35. . .0.40 although they seem to prefer larger values of v towards 
more negative values of 7; this may be explained by later onset of scaling. 
Allowing for a smooth background (in the neighborhood of the peak) like a 
constant added to the specific-heat improves the consistency of the scaling. 

Even the pseudo critical positions — although with less predictive power 
due to the uncertainty of the relation of the shift exponent A to the inverse 
critical exponent 1/V — show scaling with values of 1/ A at least roughly of 
the same magnitude. 

In general we find better scaling behavior than for the more "edgy" lat- 
tices types SH studied in [2|. Corrections to FSS due to the curvature (which 
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may be of 0(L 2 )) do not seem to be necessary. 

5 Conclusion 

After 17 years pQ of development it has now been possible to obtain a consis- 
tent picture of the scaling behavior in the pure compact U(l) lattice gauge 
theory at the confinement-Coulomb phase transition. Using a radically new 
kind of finite lattices, modern FSS methods and substantial computer power 
we have found a consistent picture of the scaling behavior of several bulk 
observables. Within the limits of numerical evidence, our analysis strongly 
suggests the existence of a non-Gaussian fixed point with v distinctly differ- 
ent from 1/2 (Gaussian value) or 1/4 (1 st order transition). Its universality 
class extends in the /3 - 7 plane along the phase transition line at negative 
7, but includes also the Wilson action, 7 = 0. 

This implies that using RG methods, one can construct a unitary contin- 
uum field theory in 4D which is neither asymptotically free nor trivial. As 
we point out in [3], this holds also for several theories related to the U(l) 
lattice theory by duality transformations. Thus rather than being an exercise 
ground for lattice QCD, the pure compact U(l) lattice gauge theory at its 
phase transition defines a sort of quantum field theory in 4D which is not 
used either in the Standard Model nor in its presently known extensions. 

The natural question is whether these novel features of the U(l) lattice 
theory might be related to the difficulties encountered in its numerical inves- 
tigation. The strongly interacting monopole loops are an obvious candidate 
for concern. The previous studies using the surface of the 5D cubic lattice [2] , 
as well as our present results, suggest that the topology of the finite lattice 
is crucial. For 7 < the two-state signal vanishes on lattices with spherelike 
topology. However, this fact alone does not yet confirm the speculations that 
the winding monopole loops are the culprits. A more detailed information 
about the field configurations, and more experience with various lattices and 
boundary conditions are required. It could be that the question is indeed 
more than technical and its pursuit might lead to a deeper understanding of 
the non-Gaussian fixed point. 
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YMP of HLRZ Jiilich and the Parallel Compute Server of KFU Graz. The 
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Appendix: Lattice details 



The lattice types SH[N] (the surface of an iV 5 hypercube) and S[N] have the 
same link connectivity structure. The lattice is build from plaquettes with 4 
links each. Not all sites have 8 links, not all links are bordering 6 plaquettes, 
and not all plaquettes are faces of exactly four 3D cubes, as it is the case for 
the hypertorus. All 3D cubes are bordering exactly four 4D cubes. 

Let us denote the total number of sites by n s , and the number of sites with 
i links by n S;i . A corresponding notation holds for the links and plaquettes. 
We have 



sites: 


n s 


= 10(N - 


l) 4 + 20(A^ - 


iy + 2 




n s ,5 


= 32 








n s ,6 


= 80(iV- 


1) - 80 






n s ,7 


= 80(iV- 


if - im{N - 


- 1) + 80 




n s ,s 


= 10(iV- 


l) 4 -60(iV- 


I) 2 






+80(iV- 


- 1) -30 




links: 


ni 


= 40(iV - 


l) 4 + 40(iV- 


I) 2 






= 80(N- 


1) 






ni, 5 


= 160(N- 


- I) 2 - 160(7V 


-1) 




ni,6 


= 40(iV - 


l) 4 - 120(iV- 


-I) 2 






+80(iV- 


-1) 




plaquettes: 


n p 


= 60(iV- 


1) 4 + 20(A^- 


I) 2 




n P,3 


= 80(iV- 


l) 2 








= 60(iV- 


1) 4 -60(A^- 


I) 2 


3D cubes: 


n 3c 


= A0(N - 


l) 4 




4D cubes: 


n 4c 


= 10(N- 


I) 4 





For example, for iV = 12 one has n s = 148832 ~ (19.6 4 ) ~ V. 

The number of plaquettes with just three 3D cubes (as well as the cor- 
responding numbers for links and sites) is suppressed relative to the leading 
terms in 0(1/N 2 ). This is typical for contributions due to curvature. We 
may say that the lattice becomes locally flat with 0(1/N 2 ). 

Contrary to the usual hypercubic torus this lattice is not self-dual. Possi- 
ble monopole loops live on the dual lattice SH', which does have a few (n p 3 ) 
plaquettes of 3 links in addition to the usual ones. 

Euler's relation for spherelike lattices (of the type discussed, i.e. without 
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further holes) is 

n s -ni + n p - n 3c + n 4c = 2 (37) 

(whereas it is zero for the torus). 

The lattice SH — in analogy to the 2D situation — may be imagined 
as an ensemble of 10 hypercubic lattices, glued together at their boundaries. 
The lattice S is constructed by a projection of SH from its center onto the 
unit sphere S 4 . 
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